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INTRODUCTION 





This Mathematics Student Book was written to enrich your study of mathematics 
by showing you how to use a computer as a problem solving tool to model mathe- 
matical concepts. The computer is particularly helpful in quickly performing the 
repetitious steps of logarithms, thus making mathematical investigations easier 
and more exciting. You will write computer programs which will help you to 
understand the major concepts involved in the study of a particular mathematical 
area. If you become more involved in investigating the laws of mathematics, these 
Books will have achieved their aim. 


To use the Lab Book for Linear Equations, you will need the following: First, 
you should have one year’s background in algebra. Second, the Lab Book assumes 
that you already know how to write a program in the BASIC programming 
language, and that you understand programming techniques for inputting data, 
performing algebraic operations, designating variables, assigning values to vari- 
ables, looping, and printing results. If you do not have this background, you will 
want to study BASIC before attempting this material. Consult the BASIC manual 
for the computer you use. Last, in order to complete the exercises in this Lab 
Book, you will need to have access to a computer for at least two hours per week. 
If more time is available, you may be able to experiment further on your own, 
either to improve your program or to investigate other areas of mathematics that 
interest you. 


Each section of this book is organized in the same way. First, the mathematical 
concepts needed to complete the exercises are reviewed. References are listed at 
the end of each section in case you want to study these concepts in greater detail. 
Next, each exercise is presented. Finally, an approach is suggested in the Problem 
Analysis and a flowchart is included to illustrate this approach. The suggested 
approach was chosen because it brings out the concepts which are being stressed, 
but the program can sometimes be written more efficiently. Once you have 
completed the exercise by following the logic in the flowchart, you are encour- 
aged to rewrite the program using more sophisticated programming techniques. 
You might also want to impose more conditions on the problem to make it more 
interesting to solve. 


There is no one “right”? way to solve a problem by programming. Experiment and 
learn as you go. You'll find you are learning something new each time, both about 
your subject matter and about using the computer to solve problems. 


Introduction 
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LIST OF SYMBOLS 





FLOW CHART SYMBOLS 


hemes Start or Stop 


ai ee dearer es Defines a process 


Se ee ee Represents computer input 


ee cet pees Represents computer output 


<> sale tun ene sr Represents a decision point 


Used to connect one part of a flow 
C3 chart continued at some other place 
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ALGEBRAIC NOTATION AND EQUIVALENT BASIC LANGUAGE SYMBOLS 


Algebraic Notation 


+ 


“or xX 


+ or / 


BASIC Notation 


/ 
SQR(x) 
ABS(y) 


INT(x) 


#or<> 


()orf[ ] 
A(1) 
A(I,J) 


RND(X) 


Meaning 
Addition 
Subtraction 
Multiplication 
Division 
Square root of x 
Absolute value of y 


Greatest integer less than or 
equal to x 


Equals 

Does not equal 

Less than 

Greater than 

Less than or equal to 

Greater than or equal to 
Replaced by 

Inclusive brackets or parentheses 
Subscripted variable 

Double subscripted variable 


Assign a random number to the 
varlable X 
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List of Symbols 


NUMBER SET DESIGNATIONS 


N - Natural number set Q@ - Rational number set 
W - Whole number set Z - Irrational number set 
I - Set of integers R - Real number set 
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LINEAR EQUATIONS 





General Form 


A linear equation defines a set of points whose graph is a straight line. The general 
form of a linear equation is 


Ax + By +C=0, 


with A40 or B40. If A=0O and B+0, the line defined by the equation is 
horizontal. If B=0 and A#0, the line is vertical. 


Linear equations are used continuously in the study of science and mathematics. 
For example, the relationship between the readings on fahrenheit and centigrade 
thermometers is defined by a linear equation, F = 9/5 C + 32° (or, in the general 
form, 5F - 9C -32° = 0). 


Slope-Intercept Form 


The general form of a linear equation can be shown to be equivalent to the 
slope-intercept form 


y=mxt+b, 


sok aU: 
where m = B and b B- 


We'll define m to be the slope of the line, determined by dividing the change in y 
(called Ay) by the change in x (Ax) as we move from one point to another on the 
line. In other words 


Ay Y2~-Yi1 
m —— rr — 
Ax X>~-Xy 





for two points (x,, y2) (x;, y;) on the line. 


Linear Equations 


(1) 


(2) 
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Let’s graph the line defined by (2), assuming that m and b are both greater than 0. 





Figure 1. Graph of y=mx+b 


The point where the line crosses the y axis is called the y-intercept, in this case 
(o,b). So you can see why we call the equation y=mx+b the slope-intercept form. 


If x, - x, = 0, then the slope is undefined and the graph of the equation is a 
vertical line. If y, - y, = 0, then the slope is zero and the graph is a horizontal 
line. 


A line doesn’t necessarily have a positive slope. Take the equation 2x + 3y -6 = 0. 
To put this in slope-intercept form, we solve for y, finding that y =- 2/3 x + 2. 
From our definition, we assume that the slope equals - 2/3. We can verify that by 
graphing the function and calculating the slope. 


° 
On alx< 





Figure 2. Graph of y = _2/3 x +2 
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Linear Equations 


From the coordinates of the points P, and P, we calculate the slope, 


Ay Yi7 ¥2 _ 4-0 2 


—— El a 


Ax xe % -3-3 3 


Point-Slope Form 
The point-slope form for a linear equation, also equivalent to the general form, is 
y-y,; =m(x-x,), (3) 
where m is the slope, (x,, y,) is a given point and (x,y) is any other point on the 
line defined by the equation. 
Two Point Form 


A linear equation can also be expressed by using any two points, (x,, y, ) and (x,, 
Y2), X; #X,, in what is known as the two-point form of the equation: 


Y2 - Yi 


yey A) (4) 
X. - Xj 


Exercise 1 — Converting the General Form to Slope-Intercept Form 
Write a computer program that will receive a linear equation of the general form, 
Ax + By + C=0, A #0 or B #0, and output the slope, the y-intercept and the 
slope-intercept form of the equation. Apply your program to the following 
equations. 

a 3x-2yt+7=0 

6b. .28x+.1ly-.6=0 

c 8x-2y=7 

dad. 8x+t+4=0 


e 2y-5=0 
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INPUT 


COEFFICIENTS A, B, C OF 
THE LINEAR EQUATION 








COMPUTE 
SLOPE, m 
y- INTERCEPT, b 






SLOPE YES 
UNDEFINED 
NO 

YES 
NO 








OUTPUT 


CORRESPONDING 
EQUATION, SLOPE 
AND y- INTERCEPT 





STOP 
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OUTPUT 


CORRESPONDING 
EQUATION, SLOPE 


AND y-INTERCEPT 








OUTPUT 
CORRESPONDING 


EQUATION, SLOPE 
AND y- INTERCEPT 
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Linear Equations 


Exercise 2 — Finding the General Form from Characteristics of the Line 


Write a computer program that will take a set of data, in any one of the forms 
discussed below and print an equation, in the general form, of a line that has 
those characteristics. Apply the program to: 


a. (2,-1) and (-3.2,7) f.  (3.27,-6) and Y-intercept=-5 
b. (0,8) and (-6,0) g. m=7/8, y-intercept= 7.3 

c.  (.8,-2) and (16,-2) h. m=85/2, y-intercept= 0 

d. (~1.5,2) and (-1.5,-8) i. m=0, y-intercept = 2 


e. (-7.5,-3) and m=-3 


Problem Analysis 


We can find an equation for a line if we have any pair of characteristics as listed 
is the cases below: 


Two points on the line 

A point on the line, slope of the line 

A point on the line, y-intercept of the line 
Slope of the line, y-intercept of the line 


me ON 


Consider the first possible case, where we are given two points on a line. Let’s 
look again at the two-point form of the linear equation. 


¥2- Yi 


Y-y, = —— (x-x,) ,X2 F X. (4) 
Xy - Xj 


We can show that equation (4) is equivalent to 
(¥; - Yo) + (X), - X,Y - Yi X2. + X,Y. =0 (d) 


which is in the general form with A= y, - y,;B=x, - x,; and C =-y,x, + 
X,y,. Thus, to find the general form of the equation from two given points, 
(x,,V, ) and (x,,y,), we need only substitute values into the expressions for A, B, 
and C above. 


In Case 2 we are given a point on the line and the slope of the line. We’ll follow 
the same procedure we used for the two-point case except our manipulations will 
involve the point-slope form, y-y, = m(x-x, ). This equation simplifies to mx - y 
-mx,+y, = 0, which is in the general form with A = m;B=-1;and C = -mx, + 
Yi- 
11 
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In Case 3 we are given a point on the line and the y-intercept of the line. This 1s 


equivalent to the case involving two points on the line, since the y-intercept 
identifies the second point, (0, b). 


Case 4 involves the slope and the y-intercept of the line. The slope-intercept form, 


y = mx + b, can be converted into the general form, -mx + y - b = 0, with A= -m, 
B=1;and C = -b. 


START 












INPUT 


SETS OF DATA, EACH SET 
1DENTIFIED AS k=1,2,3 OR 4 





k = case 


CONSIDER 1°° 


SET OF DATA 










DETERMINE 
COEFFICIENTS 
A,B,C 


OUTPUT 
By+C =0 
NO 
OUTPUT 
Axt+C =0 
NO 
OUTPUT 


Ax+By+C = 0 





DATA 
BELONG TO 
CASE k 









ALL 
DATA 
CONSIDERED 






CONSIDER NEXT 
SET OF DATA 







YES 


STOP 
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Collinearity of Points 
Three or more points are said to be collinear if they are points on the same 
straight line. We can determine if several points are collinear by finding the slope 
of the line segments joining one point to each of the other points in the set. If all 
segments have the same slope, the points are collinear. 
Exercise 3 — Determining Collinearity 
Write a computer program to determine if three or more points in the Cartesian 
plane are collinear. If the points are collinear, give the slope of the line deter- 
mined, and identify vertical and horizontal lines. 
Apply the program to each of the following sets of points: 

a. (-3,-2), (0,0) (6,4) (3,2) 

b. (-1,5) (3,5) (4.5,5) 


(2,-5) (-5,-10) (2,-5) (3,6) 


= 


d. (1,2) (-5,-10) (2,-5) (3,6) 


(2,0) (2,-3.5) (2,8) (2,- 7) 


& 


Problem Analysis 


First, compute the slopes to determine collinearity. If the slopes are all equal to 
zero, identify the line as horizontal. If the slopes are all undefined, identify the 
line as vertical. Otherwise, output the slope of the line. 


Linear Equations 
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START 







INPUT 


n, NUMBER OF POINTS IN THE SET 
COORDINATES (x;,y;) OF POINTS 
P.,.=1 TOn 












j< J+ 


OUTPUT 


“POINTS DETERMINE A 
VERTICAL LINE” 







SLOPE 
OF P; P, EQUAL 
SLOPE OF P;P.,., 


NO 






SLOPE, m, 
EQUAL ZERO 













OUTPUT 


“POINTS DETERMINE 
A LINE WITH SLOPE m” 






ALL 
SETS OF 
POINTS BEEN 
CONSIDERED 


NO 


YES 





“POINTS ARE NOT 


COLLINEAR” 








“POINTS DETERMINE 
A HORIZONTAL LINE” 
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Linear Equations 


Length and Midpoint of a Line 


Figure 3 shows an arbitrary line segment P, P,. We’ve drawn P, P; parallel to the 
x-axis and P,P, parallel to the y-axis. You can see that P; has coordinates (x,, 
y,). The length of P,P 3 is obviously x, - x,, and that of P,P; is y, - y,. 
Therefore, by the Pythagorean Theorem, the length 


P,P, = +f (x2-x)? + (yo-yi) (6) 


To determine the coordinates (x,y) of the midpoint, Py, draw lines from Py, 
through the midpoints and perpendicular to line segments P,P, and P,P. You 
can see from the figure that 

KX, +X, Vi*¥2 


5 and y = 9 











> 






P.(x2,¥2) 


Pi(x1.y¥1) ts ee aw mM P3(x>,¥1) 


Figure 3. Graph of An Arbitrary Line Segment 


Exercise 4 — Finding the Length and Midpoint of a Given Line 


Given two points (P,(x,,yV,) and P,(x,,V,) in the Cartesian plane, write a program 
that will determine the length of line segment P,P, and the coordinates of its 
midpoint. 





Apply your program to the following pairs of points. 
a. (-3,5), (2,-7) 
b. (-7,-1), (8,5) 
c.  (-5,7), (3.5,7) 


d. (-4,-6.25), (-4.4,8) 
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Problem Analysis 


Determining the length and midpoint of P,P, is a straightforward process as 
indicated by the flowchart. 












INPUT 
COORDINATES (x,,y;) OF 
POINT P,; AND{x ,y2} OF 
POINT P> 


0 
i 


LENGTH OF LINE CONNECTING 
P, AND P, 





2 
D < “VY (x2-x1) + (y2-v4)° M = x COORDINATE OF MIDPOINT 





2 


y COORDINATE OF MIDPOINT 


Intersecting Lines 


Figure 4 shows a line, P, P,, and its perpendicular bisector, L. To find the linear 
equation of L we need to know its slope and its y-intercept. 
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Figure 4. Line P, P, and its Perpendicular Bisector 


First, let’s find the slope, my,. It can be shown that the slope of L is equal to the 
negative inverse of the slope of P, P,, that is, 





1 soe 
mi, = or where m = slope of P, P.. 


If you are interested in studying the proof of this statement, consult one of the 
references given at the end of this section. At any rate, we can easily calculate the 
slope of P, P, because we are given the coordinates of two points, and that tells us 
the slope of L. 





Now, we need to determine the y-intercept, b, of L. In order to fine b, we will 
need to determine the intersection point, I(x;,y;), of L and P, P,. Remember, the 
exercise calls for the perpendicular bisector of P, P,. Therefore, I is the midpoint 
of P,P, and x; = (x, +x, )/2 and y; = (y2+y, )/2. By substituting these values into 
the slope-intercept form of the linear equation, we have 











Yj = my, xj + b. 
Since we know the values of y;, my, and xj, we now have the y-intercept, which 
we'll call by. Thus, the linear equation for L is y = myx + by, 
Exercise 5 — Equation for the Perpendicular Bisector of a Line 
Write a computer program to find the slope-intercept form of the equation of the 


line which is the perpendicular bisector of a line segment determined by two 
points P, and P.,. Apply your program to the given sets of points. 


Linear Equations 


(7) 
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a. (-7,2);(5,-6) 
b. (7,0); (7,17) 
c. (-11.5, -3); (1-3) 


Problem Analysis 


The equation for the perpendicular bisector is found by following the procedure 
discussed above and outlined in the flowchart. 


START 







INPUT 
COORDINATES (x,,y) OF 
POINT P, AND (x>,y>) OF 
POINT P5 






(x;,y;) <— MIDPOINT 
OF P,P 





YES 


YES 


OUTPUT 


y = my x + by 





STOP 
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L = PERPENDICULAR BISECTOR OF P,P, 
mp = SLOPE OF PERPENDICULAR BISECTOR 


b) = y-INTERCEPT OF BISECTOR 


OUTPUT 


CORRESPONDING 
EQUATION 


OUTPUT 


CORRESPONDING 
EQUATION 
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Another interesting problem involves finding the point of intersection for any two 
non-parallel lines. In order to do this we need to find the common solution for 
the equations of the two lines. If there is no common solution, the lines do not 
intersect. If there is more than one common solution, the equations define what 
we call co-incident lines rather than intersecting straight lines. The method of 
determining the common solution is presented in the Problem Analysis, but if you 
are interested in the details of its development, you should consult one of the 
suggested reference books. 


Exercise 6 — The Intersection of Two Lines 


Write a computer program that will find the point of intersection of two lines, 
given their equations in the general form. Apply the program to: 


a. 2x-yt3=0 
xt4yt+2=0 

c. 2xt2y=99 
xty-9=0 

b. 1/4xt+y=2 


x/2t+3y/5=14 


d. xt2y-3=0 
Axt+8y=12 
Problem Analysis 


Finding the point of intersection of two lines is a problem of finding the common 
solution of two linear equations, if such a solution exists. Lets assume the two 
lines L, and L, being considered have equations A,x +B,y+C, =Oand A,x+ 
B,y + C, =0, respectively. Further, we assume the graphs of the two equations 
are as shown below 





Figure 5. Two Intersecting Lines 


Linear Equations 
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To find the point of intersection, (xj,yj), we must solve the system below, 
i.e., find a common solution for the two equations: 


A,x+B,y+C, =0 


A,x+B,y+C, =0 


It can be shown that the system in (8) is equivalent to: 


* AB, - A,B, 
A,C, = A,C, 
* A,B, - A,B, 


Therefore, the point of intersection is 


B,C, - B,C, A,C, - A,C, 
oy 
A,B, - A,B, A,B, - A,B, 


Let us analyze what the characteristics of solutions for these equations might be. 
Let N, = B,C, - B,C, D, = A,B, - A,B, 

Ny = A.C, - A,C, Dy = A,B, —- A,B, 
Then (Nx/Dx, Ny/Dy) yields no solution if Dy and Dy equal zero and Nx # 0 
or Ny + 0. This means the equations define parallel lines. 


If Dy; Dy, Ny, and Ny all equal zero, then the system has more than one common 
solution and the equations define co-incident lines. 


Make sure you understand why there is no unique solution in the situations 
described above. 
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START 





INPUT 
A,,B,,C, A 7,B2,C, 








OUTPUT 
SOLUTION (x,y} 









A,B> = A,B, = 0 








OUTPUT 
“NO SOLUTION” 









B,C, - B.C; <9 
OR 
BiG) = Aly 0 







OUTPUT 
“MORE THAN ONE SOLUTION” 






NO 


STOP 


Exercise 7 — Finding the Center and Radius of a Circle Determined by Three 
Points 


Write a program that will consider three points in the Cartesian plane. If they are 
not collinear, find the center and radius of the circle determined. Apply the 


program to: 
a. (2,3); (7,6); (5,-4) 
b. (6,1); (-2,1); (4,8) 
c. (7,2); (-3,-5); (7,-1) 


d. (3,1); (3,-1); (3,6) 
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Problem Analysis 





Figure 6. Graph of a Circle Determined by Three Given Points 


This problem is solved by applying the techniques we have developed in the 
previous exercises. Figure 6 is a graph of the circle determined by points P,, P,, 
and P,, with some dotted lines added to help us visualize the tasks to be 
accomplished. The problem is solved by completing the following three tasks: 


(a) Find the equations of lines m and n, which are the perpendicular bisectors of 
line segments P, P, and P, P3. 








(b) Find the coordinates (h,k) of the point of intersection, 0, of lines m and n. 
This point of intersection is also the center of the circle. 


(c) Find the length, r, of segment OP, . This length is the radius of the circle. 
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START 













INPUT 


COORDINATES OF 
P,,P> AND P3 





ALL SETS OF 
POINTS BEEN 
CONSIDERED 





P, ,P2,P3 
COLLINEAR 









FIND EQUATIONS OF 
PERPENDICULAR 
BISECTORS m AND o 









(bh, k) <~ POINT OF 
INTERSECTION, 0, OF m 
AND n 


r << LENGTH OF OP, 


OUTPUT 
(h,k} AND r 


STOP 
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LINEAR SYSTEMS 


A linear system is a set of linear equations consisting of two or more equations 
with the same number of variables. The number of equations in the system is 
called the order of the system. The following system is an example of a third 
order system: 


2x+3y- z+2=0 
-x -2y+3z+3=0 
38x+ y-22-5=0 


To solve a system of equations we need to find a solution that will satisfy all 
equations in the system. The solution of the above system, if it exists, is a triplet 
(x,y,z), where x,y and ze R. 


Many problems in science, engineering, economics, etc., require solving systems of 
equations. For example, some of the systems encountered in firing guided missles 
or designing atomic reactors involve more than 10,000 variables and equations. 


In Exercise 6 we solved a second order system, using what is known as a “‘Closed 
Form” solution. By this we mean a solution obtained by a finite number of 
operations. You will remember that the form developed for a second order system 
is quite complicated, so imagine how complex a closed form for solving a 1000th 
order system would be. Even closed forms for 3rd and 4th order systems are 
unwieldy. In this section, we will learn several methods for solving higher order 
systems. 


Iterative Method 

The iterative method described below gives us an approximate solution to a linear 
system, which is adequate for most real life problems. Actually, in most problems 
the equations themselves are approximations since the coefficients are approxi- 


mate readings from instruments. 


We'll illustrate this iterative method by solving the following system: 


8x+7y=21 3 f(x) =-3/7x+3 or f(y)=-"/38y+7 


xty=2 >) g(x)=-x+2 ~~ or gly)=-yt+2 


Linear Systems 


(10) 


(11) 
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We begin by making an educated guess at the solution of the system. The actual 
solution is an ordered pair (X9,YQ). Figure 7 is a rough graph of our system. 





Figure 7. Graph of the Linear System 


The point (x,,y,) is our rough guess at the solution. In this case, we said that x, 
= 1, Substituting x, into the first equation, we come up with y, = - 3/7 (1) +3= 
18/7=2.571. Next we substitute y, into the second equation to find x, =2-18/7 
= ..571. This iterative process is continued in the table below. Notice that the 
differences in successive values of x and in successive values of y decrease as the 
approximations approach the actual solution. 


MS: eV *1.0 | -.571 |-1.245 | -1.534 | -1.657 | -1.710 | -1.733 |-1.743] -1.747 
Y, = 3/7x, 43 [2.571 | 3.245 3.534] 3.657) 3.710] 3.733] 3.743 | 3.747] 3.749 


"x4 is our initial guess, not a calculated value 









In other words, as illustrated in Figure 7, each successive point determined by the 
iterative process is closer to the solution (x9,yo). Had we continued the process 
we would have found the actual solution to be (x,7,Y,7) = (-1.75, 3.75). Our 
example thus has a rational solution, but the solution to many systems is a pair of 
irrational numbers. In that case, we would determine the amount of error 
acceptable and continue the iteration until a solution was found within the 
acceptable range. 
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For the exercises in this unit, we will accept (xy,yy) as an acceptable approxi- 
mation of the actual solution (x9,Yo9) when | Xp - Xp_-1! <e and | Yn - yn -11<e 
for some predetermined small value of e. 


In order for the iterative process to work, you must be careful to substitute the 
initial guess, x,, into the equation whose slope has the smallest absolute value. 
That is, if m,; and m, are the slopes of two linear equations (1) and (2) 
respectively, so that |m,|<|m,| then the first guess x, for xX, is substituted into 
equation (1) to find y,. You will understand why if you draw sets of intersecting 
lines with both equations having positive slope, with both equations having 
negative slopes, and with each equation having a slope of different signs. Sketch in 
the successive (xj;,yj) approximations to (X9,Y,9) when the first substitution is 
made into the line with the greatest absolute value of its slope. Do you see why it 
makes a difference into which equation you substitute your first guess? 


Exercise 8 — The Iterative Method 


Write a computer program which uses the iterative method just described to solve 
the linear systems listed below. Find solutions within a tolerance of = .0001. 


39x ty -6=0 2.1356x + 86.93y + .7394 = 0 
(a) (b) 

25x + 22y -13=0 8113x + 4.7938y - 6.8438 = 0 
Problem Analysis 


This problem is a straightforward example of the process described above for 
solving a 2nd order system of equations. 
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NO 
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T, SOLUTION TOLERANCE 
f{x) AND fly} FORMS OF FIRST EQUATION 
g(x) AND gly) FORMS OF SECOND EQUATION 





START 


INPUT 


m, <— SLOPE OF f(x) 


my < SLOPE OF g(x) 


INITIAL GUESS AT 


x, < ACTUAL VALUE 


OF Xo 










xi, — glyji_y! 


Y; < f (x;) 









¥, <— g{x;) 
|x; ~ Xj-1 I<T 
OR 
ly; —Vi-] I< T 


YES 


STOP 
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Gaussian Algorithm 


The next method for solving systems of equations involves the use of a procedure 
known as the Gaussian Algorithm. The method is based on what is referred to as 
the “elimination” or “‘addition” method in some elementary algebra books. For 
example, the following system can be solved with this method: 


(1) 2x-3y=2 
a 3x + 4y =6 


If we multiply (1) by 3 and (2) by -2 we arrive at an equivalent system 


(1) 6x-9y= 6 
(2) -6x- 8y =-12. 


Replacing (2) with the sum of (1) and (2): 


(1) 6x- 9y= 6 
(2) -17y=-6 


Multiplying (2) by - 1/17: 


(1) 6x-9y=6 
(2) y = 8/17. 


At this point we know that the second element, y, of our solution is 6/17. 
Substituting into (1) we find x = 26/17. Therefore the solution is (26/17, 6/17). 


The term “elimination” comes from the fact that in one of the equations we 
eliminate one of the variable terms. The x term was eliminated as we replaced the 
second equation with the sum of the two equations. The term ‘Addition or 
Subtraction” comes from the fact that we added two equations to get a new 
equation with which to replace equation (2). At any rate, this is a simple example 
of the use of the Gaussian Algorithm. 
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Before we can apply this algorithm to higher order systems of equations, we need 
to make a generalized statement of the procedures used above. 


The following operations can be performed to transform a system of equations 
into an equivalent system. 
a. Multiplying any equation in the system by a non-zero constant. 


b. Replacing any equation in the system by an equation which is the sum 
or difference of two equations in the system. 


c. Interchanging positions of equations within the system. 


The following third order system can be solved by applying these three 
operations: 


(1) 2x + 8y - Zz = 2 

(2) -X - 2y + 97 = 9 

(3) 3x + 3/2y - 22 = 1 
Step 1: Multiply (1) by 1/2. 

(1) x +3/ay - Voz = 1 

(2) -x - 2y + 22 = -8 

(3) 3x + 3/2y - = 1 


Step 2: Replace (2) with the sum of (1) and (2). 


(1) x +3/2y - Voz = 1 
(3) 38x + 3/dy - %” = 1 
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Step 3: Replace (3) with the sum of (3) and -3 times (1). 


(Tr). Xo 3/2y x Woes 
(2) - l/ay + 3/23 = -2 
(3) - 8y - Vaz = -2 


Step 4: Multiply (2) by -2. 


(1) x + 8/2y - Voz = 1 
(2) + y- 8 4 
(3) - 8y - Vaz = -2 


Step 5: Replace (3) with the sum of (3) and 3 times (2). 


(1) x + 3/ay ~ Woz = 1 

(2) y- 3= 4 

(3) - 19/27 = 10 
Step 6: Multiply (3) by -2/19. 

(1) x + 3/a2y - l/gzg = 4 

(2) y - 82 = 4 

(3) z = -20/19 


The equivalent system shown after Step 6 is sometimes referred to as the 
triangular form of the system. 


We now have the value of z for the desired solution. We can find the values of x 
and y by “back solving,” that is, by substituting z into (2) to find y and then 
finding x by substituting y and z into (1). 
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Back Solution: 


20 
we 3 (ea 
(2) y a6) 
76 60 
ae gee 
19 19 
_ 16 
y 19 
16 20 
3/9f — |} _1/9(-—) = 
oe (F5) (=) 1 
4 {1 _, 
oy Se 
19 19 
19 34 
2 19 | 84 
19 19 
15 
= 2? 
19 


-15 16 -20 
Therefore, the solution es, 
ererore e€ SOLUTION 1S ( 19 19 19 


The solution is not difficult to find using this algorithm for 2nd and 3rd order 
systems, but it’s obvious that the procedure will quickly become cumbersome as 
we go on to higher order systems. What we need is to find a way to use the 
computer to perform the process for us. 
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Look again at the entire process, Step 1 through Step 6. The process involved a 
series of operations performed on the numerical coefficients. If we make an array 
of coefficients of the above system, we have 


9 8-1 2 
a 2o 2 23 
g 3/2 -2 1 


which we will refer to as the coefficient matrix. The last column contains the 
constants from the right side of each equation. 


Compare the following arrays with the corresponding steps above. 


3 1 3 1 
Step 1 1 — -— 1 Step 2 1 — -— 1 
ey 2 2 oa 2 2 
1 3 
-1 -2 +2 -3 0 -— -— -2 
2 2 
3 3 
3 = -2 1 -— -~2 -1 
2 : 2 
3 1 3 1 
Step 3 1 3 “5 1 Step 4 1 ie “5 1 
0 : 2 2 ) 1 3 4 
2s ; e = 
1 1 
0 -3 -— -2 0 -3 -— ~2 
2 2 
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3 1 3 1 
Step 5 1 5 “5 1 Step 6 I 5 “5 1 
0 1 -3 4 0 1 -~3 4 
19 20 
) =~ Jf 0 1 ——— 
0 5 0 0 19 


The arrays above show that we solve an nth order system by making an n by 
(n+1) array of the variable coefficients and the constant terms. Then by applying 
the transformation principles, we try to arrive at an array in which the variable 
coefficients on the main diagonal have been changed to 1 and all elements below 
that diagonal are zero (as in the array above for Step 6). 


In general, given a system 


AyiX, + Ay2X. + Ay3x3 +°°AinXn = &, 
AyiX, + Ag2X. + Al3k3 + A,nXn = k, 
A3,X, + A32X. + A3g3x3 +°°AzgnxXn = &z 
AniX; + AnoX. + An3xX3 +°"AnnXn = Kn 


the coefficient matrix is 


Ai: Aira Ai3 «2 Ain Ai(n4t1) 
Ar, Ar2 Adz ws Aon Aj(nt1) 


A3,; Az. A33  ... Asn A3(nt1) where 


Ani1 Anz Anz... Ann An(nt:) 
34 


Ai(n+1) = K, 
Az(nt+1) = K, 
A3(nt+1) = K; 
An(nt+1) = Kn 
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We have used Ajj for the coefficients of the variables to identify which elements 
of the array are involved when one of the transformation principles is applied to a 
particular equation of the system. Also, note that the constants of the equations 
have been identified by Aj(n+,). We used xj’s to name the variables so we can 
have greater than a 26th order system if needed. 


The following array reflects the transformation of the above coefficient matrix: 


1 Fi. Ty3 
0 1 T. 3 
0 0 1 
0 0 0 


Tin Ti (nt) 


Ton T2(nt1) 


T3n 1 3(n+1) 


1 Tn(n+1 ) 


The Tjj’s are transformed elements resulting from the application of the transfor- 


mation principles. 


We will go through the steps necessary to solve a general third order system: 


AiiX, + Ajy2X 
Ay 1X + Ay 2X2 


A31X, + A32X> 


First set up the array 


Ay Ay Ay3 
Ar, Az. Ag; 


A3\ A3> A33 


+ A; 3X3 


+ A, 3X3 


+ A33X3 


Ai4 
Ada 


A34 


where 


Ky 
k, 


k 


Ai4 
Ay4 


A34 


i 


k, 


k; 
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1 
Step 1: Multiply first row by A. Aii # 0 


11 











Aj 
1 Ti2 Ty3 Tia a A. > a 
11 
Ay 
A, A>» Az 3 Axa Lis = A 
11 
Ai3 
As, A32 A33 A34 Ti3 = A 
11 
Aig 
dP A 
11 


i Tis. Tyee Lye T,, = -A,, (1) + Any = 0 
0 Too. “Toa Toa Tr. = -Ayi (Ti2) + Aor 
A3, Az. A33 Aza T,3 = -Ayi (Ti3) + Ads 

Tog = -Agi (Tia) + Aaa 


Step 3: Multiply row 1 by -A;, and add to row 3. 


1 Tye, Tye Tia T3, = -A3, (1) + Az, = O 
0 T,. T23 Toa T,, = -A;, (T,2) + As2 
0 T3. T33 T3, T3, = -A3,(T,3) + <Asj3 

T3, = -A3s, (Ti4) + Aza 
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1 
Step 4: Multiply row 2 by — 
T, 


1 Ti T 13 Ti 4 


0 T 32 T 33 T34 


5 Ts 2 


# 0 


t 
T | 


i 
T 2 


' 
Ty 3 


/ 
Ty 4 


Ty 
T > 2 





T, 2 
Th 





Ty 3 
Tyo 


T 4 
T. 2 





Step 5: Multiply row 2 by -T;, and add to row 3. 


1 T 2 Ty 3 Ty 4 


1 
Step 6: Multiply row 3 by —— 
T's 3 
1 Ti. Ti3 Tia 


0 1 T23 Too 


0 0 1 aT 


“we 


Ta 3 


T31 


T32 


T3 3 


T34 


-T 3. 
-T 32 
-T 3. 


-T3> 


J 


31 


T3 3 





T3 2 
T3 3 





T3 3 
Ts 3 
T34 
Ts 3 





+ 


(T21) 


+ 


(T22 ) 


+ 


(T2 3) 


+ 


(T3 4) 
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0 
T3, = O 
T3 3 
T3 4 
= 0 
= 0 
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In order to be able to see the “back solution” process, lets transform the array 
into the system of equations it represents. 


(1) x, + T,2X, +e “Tigxs = Tye 
(2) X2 + Ty3x3 = Toa 
(3) aa T3 4 


From (3) we see that x3 = T3, 
Back solving, we get x, = T,4 - T3 (x3) 
And x, = T,4 - T,3(%3)- Ti 2(x2) 
From the above, we can generalize the back solution for an nth order system: 


Compute successively the values of the variables beginning with xy through xj 
where the ith variable 


x} = Tint.) - Tinl¥n) - Tiqn-1)€8n-1)- °° -Tiqi+s )C¥i+1) - 


To computerize the process above, we will replace each Tjj with Ajj. As each 
transformation principle was applied above, we could have continued to call each 
new element value Ajj, but then it might not have been clear that the element 
values of the array were changing from step to step. 


Therefore 


x = Aint: )¥n+1-Ain(&n)-Ai(n-1 en- 17 777 ~Ai(its )Oi+1) » 


where 


Xn+) a 1. 
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Exercise 9 — Applying the Gaussian Algorithm 


(a) Write a computer program that will solve an nth order system using the 
Gaussian Algorithm and the back solving process described above. Apply 
your program to the systems of Exercise 8 and to the systems 


2x - y t+ 82 - 1 = 0 3x + 7y - 82 - w-2=0 

a. -x + By - 22 +0=0 4x - 2v+ 2z2- 8w-5=0 
xtily+ z2+2=0 . -2x + Qy - 52 + Qw+1=0 

4x - yt 3z+ w-38s:2=0 


Problem Analysis 


(a) The order of performing the transformation is important, and should be done 
as outlined in the flowchart below. The suggested references discuss why the 


order is important. 


INPUT 


n, ORDER OF THE SYSTEM 
ARRAY OF COEFFICIENTS, ROW BY ROW 











pen 


=") 


MULTIPLY ROW: BY 


THE RECIPROCAL OF Ai, 








MULTIPLY ROW i BY -Ax; 
AND ADD TO ROW k 





YES 
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x = Aitsy) — Ain n) — Aitne iy nad > 7° ~Aigitry) O41) 


. <a> 
YES 





(b) Adjust the computer program for part (a) above so that it will transform the 
coefficient matrix into one of the form 


I 0 
0 J] 
0 O 
0 O 
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(b) The desired matrix is obviously the coefficient matrix of the system: 


X3 


i 
nn 
tw 


Xn = Ky 


Therefore, given a system of equations, if we can transform its coefficient 
matrix to an equivalent matrix where we have the coefficients, xjj, i= 1 to n, 
all equal to 1, and all other variable coefficients equal to 0, we have the 
solution to the system. This method of solving systems of equations, which 
eliminates the back solving process, is known as the Gauss-Jordan elimi- 
nation method. 


The flowchart from part (a) must be adjusted to accomplish the following: 


1. Establish zero values for every element of each row, i, except the 
element Ajj, which should be equal to 1. 


2. Eliminate the “back solving” loop. 
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k <k+1 
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START 






INPUT 


n, ORDER OF THE SYSTEM 
COEFFICIENT MATRIX, 
ROW BY ROW 






joey 
i] 


MULTIPLY ROW, i, BY THE 
RECIPROCAL OF Aj, 





NO 


: | 


YES 

NO <> 
YES 

NO 


y 


YES 






OUTPUT 
Aint) FOR i=1 ton 


MULTIPLY ROW i BY -Ai,. 


AND ADD TO ROW k 
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Exercise 10 — Solving Linear Systems by Use of Matrices 

Write a program to solve an nth order system of linear equations by the use of 
matrices. 

Problem Analysis 

Matrix algebra is a useful tool in solving linear systems of equations. If you have 
not studied matrix algebra, you should consult one of the suggested references 


before attempting this exercise. 


To illustrate the use of matrices for this type of problem, consider the following 


system: 
2x - 3y 
3x + 2y 


By the definition of matrix equality, the above system can be considered as 
defining the equality of two 2 x 1 matrices 


(ia) (3) 


where (2x - 3y) and (3x + 2y) are single elements of the matrix. 


1 
-2 


By the definition of matrix multiplication, the matrix equation above can be 


expressed as 
2 -3 x 1 
3 2 y -2 


| : | ee 2/13 3/13 
If we multiply this equation by the inverse of ( 9 9 ) , namely | _ 3) 2) , 
we end up with 13 “13 


2/13 3/13\ (2 -3\/x\ | on a 
-3/143 2/18/\3 2/\y -3/13 2/13 / \-2 
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which simplifies to 
7 ‘\(’ ( ne 
0 1/\y SUS Ys 


By the definition of matrix multiplication and equality we can convert this to a 


system of equations 


II 

pe 

~ 
OO 


ce 


Ox + ly = -//3 


which gives us the solution, (4/ 3, ati 3), to our original system. In other words, 
we find the solution by multiplying the matrix of the right side by the inverse of 


the matrix of the coefficients. 


INPUT 











n, ORDER OF SYSTEM 

Mat A,n x n MATRIX OF 

COEFFICIENTS OF THE 
VARIABLES 

Mat K,n x 1 MATRIX OF 

CONSTANTS OF THE 

SYSTEM 


Mat S <= Inv. (A) -k 
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LINEAR PROGRAMMING 


Linear programming is a technique used to determine the linear relationship 
among a set of variables in order to obtain the optimum of a desired objective, 
given a set of restraints, expressed as linear relations, on the set of variables. We 
will illustrate this technique by solving a simple manufacturing problem in which 
the optimum relation between two variables is desired, but this technique could 
be applied to finding a desired relation among almost an unlimited number of 
variables. 


Linear programming is a very widely used application of computer technology. 
Problems demanding optimization of resources in the fields of medicine, transpor- 
tation, manufacturing, farming, nutrition and many others are solved through the 
use of linear programming. 


Exercise 11 — A Manufacturing Problem 


The owner of a manufacturing firm has a slack period in his business during which 
part of his plant, labor and capital is not utilized to the fullest. He expects the 
condition to last through the two months of August and September before his 
regular business will again utilize all his resources. 


In order to avoid laying off people and shutting down part of his plant during 
these two months, the owner has found a market for two types of Christmas toys 
which he can manufacture during this time. He can sell each type of toy in lots of 
100. The profit is $150.00 per 100 on type A toys and $250.00 per 100 on 
type B toys. A machine which makes plastic molds, will be available. For econ- 
omy, the machine must be used more than 120 hours. It can make 100 molds 
every 5 hours for toy A and 100 molds per 15 hours for toy B. The owner has 
$1300.00 capital available for this endeavor. The capital requirements are $70.00 
per 100 to manufacture toy A and $50.00 per 100 for toy B. To prepare the toys 
for market requires two types of labor skills, assembling and packaging. During 
this period, 300 hours of assembling labor will be available. Toy A requires 8.0 
hours per 100 assembly time and toy B, 20 hours per 100. There will be 400 
hours of packaging time available. They both require 20 hours packaging time per 
100 toys. How many lots of 100 of each kind of toy should be produced in order 
to yield the maximum profit? 
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Problem Analysis 


The following table summarizes the resource data we are given: 


Resource Requirements per 100 


Now we will express the linear relationships set forth by the data. In order to 
maximize profit, we find a solution to Equation (1) below where A represents the 
number of lots of toy A to be manufactured and B the number of lots of toy B. 
























Capital 


Type of Resource 
Machine time to 


make molds 


Assembly Skill 


Packaging Skill 


Total Resource 
Available 


= 


$1,300.00 


300 hours 


400 hours 










120 5 hours 










$70.00 





8 hours 






20 hours 





The other equations are derived directly from the table. 
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(1) 
(2) 
(3) 
(4) 
(9) 
(6) 
(7) 


150A + 250B = P 


oA+ 1OoB2 


120 


7T0A + 50B <1300 


20A+ 20B< 


A 


B 


8A+ 20B< 


Z 


Z 


0 


0 


300 


400 


(P, profit to be maximum) 
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B 
(3) 
30 
(5) 28 
Lin 2° 
(4) : 
(2) 





Figure 7. A Graphical Solution of Exercise 11 


Equations (2) through (5) are plotted as solid lines on the graph. The profit line, 
Equation (1), is shown as a dotted line in three positions. Writing Equation (1) in 
slope-intercept form we have 


-150 P 
= — A —- . 
B= 350° ~ 250 


At position I, the y-intercept or P/950 is 8, therefore the profit equals $2000. In 
position II, the y-intercept is about 13, which makes the profit $3250. In position 
III, the maximum position for the line to still intersect the profitable region, the 
y-intercept appears to be 17, which makes the profit $4250. The point M seems 
to be the point with integral coordinates nearest the highest intersection point of 
the profit line. The coordinates of M are (9,11), which means 9 lots of 100 of 
type A toys and 11 lots of 100 of type B toys will yield the greatest profit. The 
actual profit should be 9(150) + 11(250) = $4100. We will write a computer 
program to determine the accuracy of our graph. 
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In writing the program we know from Equations (6) and (7) that the minimum 
value of A or B is zero. From Equation (3) we can determine that the maximum 
values are A = 18 and B = 26. We need to consider only integer values of A and B, 
because one constraint was that the toys be manufactured in lots of 100. 


30 


START 












STORE AANDB 
AND CORRESPONDING 
PROFIT 







CONSTRAINTS 
{2) THROUGH (5) 
SATISFIED 


OUTPUT 
MAXIMUM PROFIT 
AND CORRESPONDING 
AANDB 
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NOTES 


